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Abstract 



In this work, the properties of the p meson at rest in cold symmetric nuclear matter 
are studied. We make use of a chiral unitary approach to pion-pion scattering in the 
vector-isovector channel, calculated from the lowest order Chiral Perturbation Theory 
{xPT) lagrangian including explicit resonance fields. Low energy chiral constraints are 
' considered by matching our expressions to those of one loop xPT. To account for the 

medium corrections, the p couples to tttt pairs which are properly renormalized in the 
nuclear medium, accounting for both p — h and A — /i excitations. The terms where the 
p couples directly to the hadrons in the p — /i or A — ft, excitations are also accounted 
Q ' for. In addition, the p is also allowed to couple to iV*(1520) — h components. 

C ■ PACS: 14.40.-n: 21.65.+f; 12.40.Vv 

Keywords: Rho meson; Medium modification. 
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X. 

c3 ■ 1 Introduction 

One of the major goals of nuclear physics is the understanding of the properties of hadrons 
at the high baryonic densities that occur in nuclei or neutron stars. However, it is difficult 
to extract information on the mass or width of an in-medium meson from experimental 
data because in most of the decay channels the final particles undergo a strong distortion 
before they get out from the dense matter and can reach the detector. 

Electromagnetic decays offer, at least in principle, a cleaner probe of the high density 
regions. Neutral vector mesons, like , are specially interesting due to its large leptonic 
width. Any changes in the mass or width of these mesons should be reflected in the invariant 
mass distributions of the leptonic decay products. 

Although there are some problems in the interpretation of data, due to the low statistics 
and the existence of many additional sources, dilepton spectra, such as those measured at 
CERN lH, ^, ^ and at lower energies at Bevalac [Q, ^, may indicate either a lowering of 
the p mass or a large broadening of its width. In the near future, experiments at GSI 
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(HADES Collaboration [^,0]) could clarify the situation by providing better statistics and 
mass resolution. 

Many theoretical approaches have been pursued to analyse the rho meson properties in 
the medium. A good review of the current situation can be found in ref. Q. Here, we 
only present a brief description of the main lines of work. Much of the current interest 
was stimulated by Brown and Rho who by using scale invariance arguments, obtained 
an approximate in-medium scaling law predicting the p mass to decrease as a function of 
density. Also Hatsuda and Lee |10| found a linear decrease of the masses as a function of 
density in a QCD sum rules calculation, although recent works |l^ cast some doubts 

on these conclusions and find that predictions from QCD sum rules are also consistent 
with larger p masses if the width is broad enough. Several groups have investigated the p 
selfenergy in the nuclear medium studying the decay into two pions, the coupling to nucleons 
to excite a baryonic resonance, or both. The main effects on the two pion decay channel 
are produced by the coupling to A-hole excitations [0, 15, 16, 17, |l^, Whereas the 



p mass is barely changed, quite a large broadening, which provides a considerable strength 
at low masses, is obtained in these calculations. Another strong source of selfenergy is the 



excitation of baryonic resonances, mainly the A^*(1520) considered in refs. |21, 22, 23| 



which further broadens the p meson. Finally, we could mention the systematic coupled 
channel approach to meson-nucleon scattering of ref. |2^ that automatically provides the 
first order, in a density expansion, of the p selfenergy, finding also a large spreading of the 
strength to states at low energy. 

In this paper we investigate the rho meson properties in cold nuclear matter in a non- 
perturbative coupled channels chiral model. This approach combines constraints from chi- 
ral symmetry breaking and unitarity, and has proved very successful in describing mesonic 
properties in vacuum at low and intermediate energies [^]. The effects of the nuclear 
medium on the mesons have already been analyzed in this framework, for the scalar isoscalar 



(I = J = 0), "a" channel |27, 28 1. There have been several attempts to translate the nu- 
clear matter results into observable magnitudes, which can be compared with experimental 
results^, However, as mentioned before, it is difficult to disentangle the meson prop- 
erties at high baryonic densities due to the strong distortion of the decay products. The 
rho meson studied here is better suited for such a study because of its large leptonic width. 

In section 2 we give a review of the model for meson-meson scattering in vacuum, 
presenting results for the p channel. Nuclear medium corrections are studied in section 3, 
with special attention to requests of gauge invariance. In section 4 we describe how the 
coupling to A^*(1520) — h components can be incorporated in the calculation. Our results 
are presented in section 5. Finally we summarize in section 6. 



2 Meson-meson scattering in a chiral unitary approach 

In this work we study the p propagation properties by obtaining the tttt — > vrvr scattering 
amplitude in the (/, J) = (1, 1) channel. We first give a brief description of the model for 
meson-meson scattering in vacuum, and then discuss modifications arising in the presence 
of cold symmetric nuclear matter. The model for meson-meson scattering in vacuum is 
explained in detail in [32|. Following that work, we stay in the frame of a coupled channel 
chiral unitary approach starting from the tree level graphs from the lowest order 
lagrangians including explicit resonance fields |34]. The model successfully describes 



vrvr P-wave phase shifts and vr, K electromagnetic vector form factors up to ^/s < 1.2 GeV. 
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Low energy chiral constraints are satisfied by a matching of expressions with one-loop X-f ^• 
Unitarity is imposed following the N/D method. 

We start from the (/ = 1) vrvr, KK states in the isospin basis, using the unitarity 
normalization [31|: 

\KK) = ■^\K+K- -K^K^). (1) 
V2 

Tree level amplitudes are collected in a 2 x 2 matrix whose elements are 



i^ii(s) 

Ki2{s) 
K2i{s) 

K22{s) 



3 P 
Kuis) 

3/2 



2Gl 



1 + 



2G2 



P Ml-s_ 



1 



2Gl 



P Mj,-s 



(2) 



with the labels 1 for KK and 2 for vrvr states. In equation (^) Gy is the strength of the 
pseudoscalar- vector resonance vertex, / the pion decay constant in the chiral limit, s the 

squared invariant mass, Mp the bare mass of the p meson and pi = s/'l — mf. In any of 

the amplitudes of eq. (^) the first term comes from the 0{p^) chiral lagrangian, while the 
second term is the contribution of the resonance lagrangian [^] . 

The approach to meson-meson scattering in the isovector channel commented so far is 
essentially equivalent to using gauge vector fields for the p meson due to VMD. We will use 
the gauge properties of the p mesons in order to evaluate vertex corrections which appear 
when self energy insertions are also considered. 

The final expression of the T matrix is obtained by unitarizing the tree level scattering 
amplitudes in eq. @. To this end we follow the N/D method, which was adapted to the 
context of chiral theory in ref. [51] . We get 



T{s) = [I + K{s)-G{s)r'-K{s), 



(3) 



where G{s) is a diagonal matrix given by the loop integral of two meson propagators. In 
dimensional regularization its diagonal elements are given by 



GFis) 



1 



16 7r2 



-2 + di+ Ui 



log 



Ti{s) + 1 



(Ti[S) 



(4) 



where the subindex i refers to the corresponding two meson state and ai{s) = y 1 ~ ^nif/s 
with rrii the mass of the particles in the state i. 

In Fig. |l| we show a diagrammatic representation of the three matrices involved, i^(s), 
G{s) and T{s). The last one represents the resummation performed in eq. (^) with the 
unitarization procedure (N/D method). 

The di constants in eq. (§) are chosen to obey the low energy chiral constraints [^ ], 
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K = 'v: + : : = ;/////// : 

(a) ^ (b) 

G = 



T = ]///////[ + ^^ZZZZZZA W///// + 



+ yzzzuA y/////A y////// + 



Figure 1: Diagrammatic representation of the matrices described in the text, (a) and (b) correspond to 
the amphtudes derived from the 0(p^) and resonance chiral lagrangians respectively. The double solid line 
filled with short oblique dashes represents the effective scattering amplitude. Dashed lines are tt, K mesons. 



-2 mi / ml. 1 , ml- 1 

dK = — ^— log ^ + - log ^ + - 

— m^ \ z /x^ z 

m% / mi 1. mi 1\ , , 

m|^ — m^ \ /i^ 2 /x^ 2 / 

and they are obtained by matching the expressions of the form factors calculated in this 
approach with those of one loop xPT. In eq. (|) ^ = 770 MeV and we have substituted 
index values i = 1, 2 by their actual meaning. 

The regularization can also be carried out following a cut-off scheme. The loop functions 
with a cut-off in the three-momentum of the particles in the loop can be found in the 



appendix of ref. |26|, and up to order -^kx they read 

(6) 



16 vr^ 



-2 log — Vai{s) log 



mi ai{s) - 

where g™"^ is the cut-off mentioned above. By comparing the expressions of the Gi{s) 
functions in both schemes we can get the equivalent 5™"^ in order to keep the low energy 
chiral information of the di constants, and therefore one can use equivalently any of the two 
procedures. 

In Fig. |2| we show the T22 matrix element vs invariant mass. To check how important 
is the inclusion of the KK channel we have performed the calculation in both coupled and 
decoupled cases. By decoupled case we mean turning off the KK in: interaction, what is 
easily done setting K12 = in eqs. (§) and (^). The conclusion is that kaon loops produce 
minimum changes in the results and therefore they can be ignored and we can work in a 
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decoupled description in what concerns the vrvr tttt scattering amphtude. In the next 
sections we shall develop the model for scattering in nuclear matter. Although we formally 
work in the general scenario of the coupled channel approach, we perform our calculations 
in the decoupled case. 



30 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 



Mj (GeV) 



Figure 2: Real and imaginary parts of the tttt — > nn scattering amplitude {T22 matrix element). The solid 
line curve is obtained with the coupled channel formalism including kaons, while the dashed curve is the 
result of the decoupled calculation. 

It is interesting to establish connection with other approaches where tadpole terms are 
explicitly kept in the lagrangian. Indeed in [^] an explicit term appears in the lagrangian 
when the minimal coupling is used to introduce the gauge vector fields. A contact term, 

>C(^*)=9>,P^^+vr_, (7) 

is obtained where gp is connected to our coupling via 

a - ^^^^ (R) 

With this lagrangian one obtains a contribution to the p selfenergy which is expressed 
in terms of the tadpole diagram in Fig. ^ and is given by 
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Figure 3: Tadpole diagram contributing to the p meson selfenergy, coming from the pp-K-K vertex. Wavy 
lines represent the p meson in the gauge vector field formalism and dashed lines are pions. 



'i'ttad 



(2vr)' 



(9) 



where DQ{q) is the free pion propagator. When we work with this formahsm we simulta- 
neously release the on shell condition in the pmr vertex in the evaluation of the selfenergy 
diagram represented in Fig. [E^, 35 1. Since we are interested in the p at rest, we need 




Figure 4: Two pion diagram contributing to the p meson selfenergy. 

only to evaluate the spatial components of H^'^ and, thus, using the vector form of the pvrvr 
coupling of 1 18 1, 



we obtain 



:iq'qWoiP-q)Doiq)-igl 



(2vr)' 



(2vr)' 



r2g'Wo{q). 



By using dimensional regularization it is easy to prove that 



(27r)4 (g2 + s + iey 



d'^q 2g' 



(2vr) 



4 „2 



+ s + ie 



(10) 
(11) 
(12) 



6 



Then eq. (|TT| ) can be written as 



Sq f dq^ 



(27r)2 



27r 



f[Do{P - q)D^{q) 



d 



drrii 



Do{q)]. 



The q^ integration is now easy to perform and we obtain 



9p 



(27r)3 



+ 



uj{q) (PO)^ - iuj{q)^ + i€ cj(g)4a;(g)2 



(13) 



(14) 



ml + q"^. 



where u!{q)'^ 

Now, if we separate in the first term of eq. ( [T^ ) the on shell and off shell contribution of 



the factor q'^ in the numerator, which comes from the pm: vertices, we have q' 



^2 



+ 



with q ; 



(^) 



-mz 



Since we have (P 



0\2 



-^u{qf 



^{q Qon) then the contribution 



of the tadpole plus the off shell piece of the two pion loop diagram is given by 

^2 r ^3 I 



(27r)3 uj{qY 



(15) 



which is a logarithmically divergent piece independent of P^ and hence can be absorbed in 
a renormalization of the p bare mass. 

The former exercise proves that the formalism keeping tadpoles and full off shell depen- 
dence of the pTTTT vertex is equivalent to the one we have showed above where only the on 
shell part of the pTrn vertex is kept and no tadpole is included. 

In the medium, however, the pion propagator in the tadpole term will change. Hence, 
in order to stick to the gauge invariance of the vector field formalism we shall keep this 
tadpole term |12, 17, 18]. 



3 (/, J) = (1, 1) TTTT scattering in the nuclear medium 

We now address the calculation in the presence of nuclear matter. As already mentioned, the 
p meson couples to intermediate two-meson states. Medium corrections will be incorporated 
in the selfenergies of the mesons involved and the related vertex corrections as well as from 
direct coupling to the hadrons in particle-hole {p — h) and delta-hole (A — h) excitations. 
We briefly discuss below the input used for the pion selfenergy. 



3.1 Pion selfenergy in dense nuclear matter 

The pion selfenergy in the nuclear medium originates from p — h and A — h excitations, 
corresponding to diagrams in Fig. |5|. We describe the selfenergy as usual in terms of 
the ordinary Lindhard function, which automatically accounts for forward and backward 
propagating bubbles, hence including both diagrams of Fig. ^. Short range correlations are 
also considered with the Landau-Migdal parameter g' set to 0.7. The final expression reads 
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where here q is the four-momentum of the pion, p stands for the nuclear density and 
C = (-^^j^)^- The Lindhard function, U = Un + C^A, contains both p — h and A — h 
contributions. Formulae for Un and J7a with the normalization required here can be found 
in the appendix of ref. [^]. We use a monopole form factor f{q'^) = fj^j^^i for the ttNN 
and ttNA vertices with the cut-off parameter set to A = 1 GeV. In the case of the ttNA 
vertex we also include a recoil factor / ^/sA to account for the fact that the momentum 
q in the ttNA vertex should be in A CM frame. We will check the sensitivity of our results 
to variations of both the A and g' parameters. 




Figure 5: Direct and crossed p — h, A — h excitation terms giving rise to the pion selfenergy in the nuclear 
medium, expressed in terms of the Lindhard function as described in the text. 



3.2 Gauge invariant contact terms 

In our microscopic model of vrvr scattering in nuclear matter all graphs in Fig. ^ are con- 
sidered. In addition to the single p — h excitation (diagram ^) we have to account for the 
medium modifications of the /?-meson-meson vertex via the /o-meson-baryon contact term, 
as requested by the gauge invariance of the theory ||, [l^]. Hence we need to evaluate 



also the diagrams ^3 to |6|d. This contact term can be constructed in a similar way as it 
was done in [^] for the 4> decay in the nuclear medium. Let us study the set of diagrams 
depicted in Fig. |^, for the case of an incoming vr^. The amplitude corresponding to the 
graph l^a in the gauge vector formulation of the p is readily evaluated and its contribution 
is 

it = z^^e.iQ + 2Kr——^^ V2i^±^)a{K + Q), (17) 



where K and Q are four-momenta of the pion and the p meson respectively, a is the spin-^ 
operator and is the polarization vector of the p meson. The contact term, represented 
by Fig. must have the structure 

- itcont = e^^^. (18) 

and gauge invariance requires that if we replace by then the sum of all the terms in 
Fig. ^ vanishes. To find out the strength of the coupling we shall take the limit — > 
because in that case diagrams Is and |d give a null contribution and we have to consider 
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/ / / / + lower line 

(a) 



IZZZr "TZZZ TZZl ^/ZZZ + lower line 

(b) (c) 



////y ^ \ S > //// + lower line 

(d) 

Figure 6: Medium correction graphs: single solid lines are reserved for particle-hole excitations. 

only the sum of the first two so far commented. Then we find the following system of 
equations: 

-^r2{^)aQ-BQ = ^. (19) 
After solving for i?^, the amplitude corresponding to diagram 03 is 



and with this information the whole set of diagrams in Fig. can be evaluated. Note that 



itcont = -^V2{^^)ae, (20) 



because of the p coupling to the pions, eq. (20) has a minus sign for an incoming vr 



3.3 Technical implementation of vertex corrections 

In the previous section we have commented that gauge invariance requires the presence of 
the p-meson-baryon contact term, which results in the need of including diagrams b, c and 
a m Fig. |. Now we pay attention on how these vertex corrections can be incorporated in 
the unitary treatment of the T matrix described in section 2. To this end let us consider 
in detail the set of diagrams ^a, b and c ^. The momentum labelling is shown in Fig. ^ 
for the first diagram (the pion in the upper line will have four-momentum q and the one in 
the lower line {P — q)). The amplitudes originating from these graphs in the gauge vector 
formulation of the p are evaluated using Feynman rules and the results are 



Aa) = I ^V,UqfiDo{q)iDo{P 



2/ 



rU{q,p) 



iDoiq) 



^We restrict here the discussion to the case of pions. The argument for kaons would be analogous. 
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fc^ (d) 



Figure 7: Gauge invariant set of diagrams used to calculate the p-meson-baryon contact term as explained 
in the text. The wavy lines represent here the vector meson, the dashed ones pions and the solid lines are 
reserved for nucleons. 



A 



^ 2/ ' 



'Uiq,p) 



iDo{q)iDo{P 



A 



A 



where V^TrTr stands for the pTrvr vertex, 



P 



2eq 



(21) 



(22) 



and DQ{q) is the bare propagator of a pion with four-momentum q. The subindices (a), (b) 
and (c) refer to each diagram in Fig. |. In A( we have exphcitly substituted the expression 
of the pion selfenergy due to a. p — h or A — h excitation, which would read 



-^U{q,p) = -^{^^ffU{Q,p)■ 
It is obvious that the sum of all three diagrams, 



(23) 



^(a) + ^(fe) + ^(c) 

d^q 



Vp^^{qyiDo{q)iDo{P - q) 



2/ 



U{q,p) 



(q^Doiq) + 1) 



(27r)4 ''"^ 

can be cast in the form of A(^a) if we substitute the q^ appearing in eq. ( p3[ ) by 

q^^q^ + D,\q). 



(24) 



(25) 



We still have to consider the last diagram in Fig. ^, It is easy to find, in the same way 
as before, that the amplitude can be cast in the structure of Af^^) niaking the substitution 



3 VM 



(26) 
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"\ Do 



/q ^ \ 

\ t 



Figure 8: Momentum labels for the in-medium correction diagrams. Do indicates a pion propagator and 
n refers to the pion selfenergy due to p — ft excitation. 



in eq. (^31). Hence, we find that vertex corrections can be automatically introduced in the 
gauge vector formulation of the problem by performing the following substitution 



f - HU<1) - f + D,\q) + (27) 



in the pion selfenergy described in eq. (pi3D, hence working with a kind of an 'extended' 
pion selfenergy which takes care of the medium modifications of the pirn vertex, namely 

(D+EYulq p) 

The calculation of the tttt scattering matrix in the medium follows the calculation of the 
G{s) function for pions. As it was mentioned in section 2 we keep full off shell dependence 
of the /JTTvr vertex. The medium effects are introduced by a subtraction of the p meson 
selfenergy due to pion loops in free space from the same quantity calculated in nuclear 
matter. The following substitution is performed in eq. (^): 

G{S) ^ G{S) + ^[/w(s) - Ifreeis)], (29) 
9 on 

where 

Imed{s) = ij -0^fD{q)D{P-q) 

■4 

(2^)4 

and D(q) is the in-medium pion propagator. 



Ifreeis) = > I l^U?' Do{t,)Do(P - q) (30) 
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To include the vertex corrections discussed above, the pion selfenergy 11^ ((7,/?) in the latter 
equation must be replaced by the one in eq. (p8[). 

Each of the integrals in eq. (^) are quadratically divergent. The subtraction Imed—I free 
cancels these quadratic divergences and we are left with a logarithmic divergence which in 
principle has to be regularized with a cut off in the pion loop momentum. However, notice 
that D{q)D{q) — DQ{q)DQ{q) is proportional to Il^^^{q, p). Therefore, since all the pieces that 
enter Il^^*{q, p) are multiplied by a squared form factor f{q'^)'^ the subtraction is convergent. 

3.4 Rho meson propagator from the T matrix 

For the sake of simplicity, and having in mind the results of section 2, we turn off the 
interaction between tttt and KK systems and work in a decoupled channel approach {K12 = 
in eqs. (I) and (D). Then the vrvr TTTT T matrix element takes the simple form 

Let us first assume that — ^ = 1, as it is indeed the case if VMD holds exactly. Then the 
K22 matrix element reads 

2 1 

K22 = -3 >%_M2' 

where one can identify the p-wave structure of the vrvr tttt scattering amplitude mediated 
by the exchange of a p meson. Hence the factor (1 + K22G)^^ in eq. ( |3^ ) is accounting for a 
selfenergy of the p meson arising from the coupling with vrvr pairs (to switch on the presence 
of the nuclear medium one has just to perform the substitution of eq. ( ]29| ) as described in 
section 3 with the modified pion selfenergy of eq. (p8[)). We can therefore write eq. (^) as 

2p^M2 1 

-t22 — 



3 p s-M^-n;-'"°p 



n^^/oop ^ _K,,G{s - M^^) = -^^G, (34) 

jTT — loOp ■ 



where Hp ^ is the selfenergy commented above. 

2q2 2G^ 

Now we consider the case -j^ 7^ 1. Let a = -j^ and proceed as above. The K22 matrix 

the " 



element can be easily cast in the form 

■3 /2 -^-'s-Mj 



2 plM^ 1 
K22 = -/-^B{s)-—, (35) 



where B{s) = 1 + (a — Eq. (|35|) reduces to eq. (P3) if a = 1, because in that case 

B = 1 independently of s. For the values of Gy and / of we find that a takes the value 



~ 0.7. For instance, at s = we have B{Mp) ~ 0.7, therefore a 30 % of deviation from 
VMD is introduced. The T matrix element for vrvr — > vrvr scattering now reads 



t22 



nr'""'' = \^B{s)G = ~g'\plG, (36) 



12 



where = ^^j^B(s) defines a coupling which will replace gp in our formulae, and Dp is 
the p meson propagator. 

In the following sections we shall discuss about the tadpole term in the nuclear medium 
as well as new contributions due to the direct coupling of the p meson to the baryonic lines 
in particle and delta-hole excitations. These new contributions will manifest as additional 
pieces to the whole p meson selfenergy 11. The way of introducing these new pieces in 
the T matrix and the p propagator follows the same scheme as it has been done with the 
contribution from two pion loop selfenergy plus its vertex corrections, 11^"'°°^. 

3.5 Tadpole term in nuclear matter 

We have seen that keeping the full off shell dependence of the pmr vertex and explicitly 
considering the free tadpole term in the calculation of the p meson selfenergy is equivalent 
to calculating with the on shell part of the pirn vertex and describing the pion loop con- 
tributions with the G{s) functions defined in section 2. This freedom allows us to include 
the tadpole diagram in nuclear matter by performing a subtraction between the selfenergies 
corresponding to the tadpole term in the medium and in free space. In the gauge vector 
field approach one has 



zU%' = -2glgp, / -^l)o(g). (37) 



(27r)^ 

The p selfenergy, 11, which appears in the p propagator in the former section can be related 
to Hpy by 

n;.. = |-<7M. + ^^|n, (38) 

and restricting ourselves to the space components of /i, v we have for the tadpole piece 

n*'^'^ = 2igl I ^,Do{q). (39) 

We define the subtracted in medium tadpole contribution to the p meson selfenergy as 

&^\p) = 2igl I ^,[D{q) - Do{q)]. (40) 

Given the analytical structure of D{q) and DQ^q) by performing a Wick rotation one 
can see that n*'^'^ is real, and eq. (^) can be rewritten as 

n*"'^(p) = -2gl I ^[ImDiq) - ImD,{q)]. (41) 

Taking into account that both D^^q) and D{q) are even functions of q^, the integration in 
q^ can be replaced by twice the integration from to oo. This integration can be performed 
analytically for the second term of eq. (^) and we obtain 



4w(g) 



(42) 



Eq. (42) is free from quadratic divergences. Actually the form factors present in the pion 
selfenergy make it convergent and a cut off in the pion loop momentum is not required. 
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3.6 Other medium corrections 



As it has been discussed in subsections 3.2 and 3.3, gauge invariance provides us with pirNN, 
pnNA contact vertices which give rise to the diagrams depicted in Fig. ^. In addition to 
those there are other medium corrections that arise within the gauge vector field formahsm 
and generate interactions via the minimal coupling scheme (the former contact vertices can 
also be obtained in this way). Let us summarize the interaction lagrangians of ref. |p!8| 
involving pvrvr, pp-mr, ttNN, ttNA, pNN, pnNN, pAA and pirNA vertices: 

^np = Ugpp^{n$d''$+d''$Ts$)-^glppp''Ts$T3$ 
C^NA = -—ipf^ippd^'^ + h.c. 

"T-vr 

CpNN = -y^M'^ 
C-pttNN = igp i!-f^ ^TTpTsfl) 

Cp^NA = -igp—i^f^i^^p^'Ts^ + h.c. (43) 

From the list above, lagrangians generating direct couplings of a /? meson to either a 
nucleon or a delta were not considered in subsection 3.2. As we will discuss below, these 
terms give rise to a set of diagrams correcting both the p meson selfenergy from two pion 
loop and tadpole terms. 

Before proceeding with the discussion of the new contributions, it is interesting to com- 
ment on the relative sign of the pmr coupling constant in the language of chiral lagrangians 
and in the gauge vector field formalism. From eq. (|4^), the piriT vertex takes the form 

- itpTTTT = -igpepip - p'Y, (44) 
with the momentum labels corresponding to Fig. However, in our model we would write 



.MpGv / /N/, 
itp^^ = i—y^e^{p - p)^', (45) 



and this establishes the correspondence 



to connect both schemes. As an example, following the lagrangians listed above, eq. ( P0| 
reads 



itcont = -gp—V2ae, (47) 



which reconfirms the relationship of eq. (|4 
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p' 



Figure 9: pn-K vertex. Lines as in Fig. 




Figure 10: Vertex correction involving pNN, pAA vertices. 
3.6.1 Extra two-pion loop correction to the p selfenergy 

The first additional medium correction that we consider here is the one shown in Fig. |l^, 
with the p meson attached to a particle line. As we shall show below, this contribution 
can also be cast in the form of a vertex correction to be introduced in the 'extended' pion 
selfenergy of eq. (p8|). Following the Feynman rules for the vertices involved (after a non- 
relativistic reduction) the p meson selfenergy contribution from this diagram (N — h case, 
po = ^) reads 

ipw(9» - 9l - UMa'', q) + UN{-q\ q] - UNi-q" + VS,q}\, (48) 

where the Un functions account for the direct contribution of a nucleon-hole excitation (Fig. 
^, left) and sum over spin but not over isospin. They satisfy the relation 

2[UN{q'',q)+UN{-q'',q)] = UN{q\q). (49) 
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In the case of a p meson coupling to a A line one obtains a similar expression, with 
different spin/isospin factors and a scaling factor to the vrvrA^ coupling, if a constant A 
width is considered. The reason is that in the nucleon case one is allowed to express the 
product of two particle propagators as the difference of such propagators times s~2 . Then 
the resulting expression of I!*-' can be written in terms of the Un functions. However, 
keeping the A width is important and, hence, we shall keep its energy dependence and use 
another approximation to the integral of the two A and one hole propagators appearing in 
the expression of . The calculation is the done by making a Fermi average over the two 
A propagators. 

Given the structure of this correction it is convenient to incorporate it by means of cor- 
rections to the pion selfenergy in the two pion propagators of eq. (^) . This is accomplished 
by means of the substitution 

fU{q'>,q) ^ fUiq'>,q) + {aN + aA)D^Hq^q), (50) 



with 



aN = A Uwiq^ -V^,q)- ^Jv(g°, q) + UN{-q°, q) - UN{-q° + V~s, q)] 



72 in / ^. \ 2 



OA 



G'A(-g°,g1GA(-<z° + V^,9l} 



(51) 



and 



GA(g°,g) 



^-Ma + |rA(sA) 



SA = + {qy - f + 2q\lMl + (52) 



Expressions for the energy dependent A width can be found in [37|.The substitution of eq. 
(|50|) already accounts for the two possible time orderings in Fig. 10: incoming p meson 
coupled to the baryon loop, outgoing p meson closing the two pion loop, and vice versa. 
Notice that in eq. (|5^) we have not used the common non-relativistic approximation of 
since it is not appropriate for this case because high momenta contribute to the loop 
integrals. The same argument is of application for the corrections discussed in the next 
section. 



3.6.2 Tadpole like corrections 

We study now a set of diagrams containing pNN, pAA vertices which are generated from 
those in the former section by contracting one pion line, in analogy to the way the vertex 
corrections in Figs, ^b, ^c where generated from Fig. |^. These terms are related to the 



diagrams labelled as ppmr vertex correction in [|1^]. The diagrams are shown in Fig. 11 
although the set includes also the partners with an outgoing p meson attached to a nucleon 
(or A) line. 

The selfenergy corresponding to the diagrams in Fig. 11 and their partners for the 
nucleon case reads 

-'^tNN = lj^^[{^yS'' I ^,q'Do{q)[UNiq'-V-s,ql-UNiq' + V-s,ql], (53) 
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Figure 11: Additional medium corrections. The pion line represent either a tv^ 
with the outgoing p attached to the nucleon (or A) line are not displayed. 



or a TT . The partners 



where the factor four indicates that the contributions of protons or neutrons over the Fermi 
sea are identical, and the same happens for the diagrams with the incoming and outgoing 
p mesons interchanged. 

For the case of the A isobar, the number of diagrams involved are twice the number for 
the nucleon case, since now the A can be excited to its four charge states. As it happened 
in the previous subsection the result cannot be expressed in terms of modified Lindhard 
functions. Instead, we use the same approximation as before for the integration of the 
fermionic loop and we get 



3 Ma 9 \m^J \fNj J (27r)4 



^\ Ga(9° + \/^,glGA((Z°,gl +GA((?°,g)GA((?° 



(54) 



3.6.3 Hartree like nucleon tadpole term 



In ref. |18| we find no other medium corrections to the /jvrvr and ppmr vertex functions, since 
the p coupling to baryons is treated in a non-relativistic approach up to order (I/Mn)^. 
At next order an additional four point ppNN vertex enters the calculation. This can be 



seen in ref. |17| where this vertex is generated by minimal substitution in a non-relativistic 
treatment of the pNN lagrangian. 

The same result can be obtained by calculating the pN pN amplitude in the rela- 
tivistic theory and then performing the non-relativistic reduction up to order M^^. This 
is schematically shown in Fig. In this non-relativistic reduction we obtain the nucleon 
direct and crossed pole terms of Fig. 12b with non-relativistic propagators and vertices. 



plus a contact term of order . 

If the nucleon lines of this last diagram are closed one obtains a nucleon tadpole dia- 
gram (Fig. [l3|b) which contributes as an additional linear-in-density piece of the p meson 
selfenergy |17|. The Feynman rule for the ppNN vertex (Fig. p!3|a) can be derived from the 
pN pN amplitude and it reads 



it 



ppNN 



1 1 



P4Mn 



(55) 



and the nucleon tadpole term is easily evaluated, 

1 f d^p 



'■'■'■ppNN 



-'"''jr. 



1 



n{p) 



(27r)4 



d^p 
(2^ 



n{p) = -iS^'gl 



£n{p} - 
1 1 



le 



4 Mat' 



(56) 
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reduction 
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% % 
(a) (b) 

Figure 13: (a) ppNN vertex, (b) Nucleon tadpole term. 

The result is an energy independent function proportional to the nuclear matter density. 
Its inclusion as a selfenergy term in the p meson propagator follows the same procedure as 
already discussed for other contributions. 

Finally, some other medium corrections involving pBB, pirBB' and ppBB vertices can 
be considered but they either vanish or are very small and we discuss them briefly in the 
appendices. 



4 A^*(1520) — h contribution to the p meson selfenergy 

Next we intend to modify our formalism so that the p meson is able to excite baryonic 
resonances via a pNN* coupling. From now on we will restrict the discussion to the case 
of Ar*(1520) I{JP) = although the same arguments considered here may be used to 

introduce other resonances with different quantum numbers 0. 

The correction under discussion (excitation of a N* — h pair) will manifest as an extra 
selfenergy term in the p propagator. The basic vertex involved in this effect is shown in 
Fig. p^a, and the lagrangian describing the interaction reads 

J^N'Np = -gN*Np^NSi4>iT'^N'- + h.C, (57) 

^Note that since we stay in the center of mass frame of the p meson, the contribution of resonances which 
couple to the p in P-wave vanishes. The coupling to A''* (1520) is S-wave and therefore its contribution to 
the p selfenergy survives even at zero p three-momentum. 
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where "^n, ^n*, 'Pi ^-re the N, A^*(1520), p fields, Si is the ^ — > | spin transition operator, 
r is the isospin-^ operator and gN*Np stands for the pNN* coupling constant that we take 
from [38| to be gN*Np = 7.73/\/3- Feynman rules from the lagrangian in eq. ( [5^ ) allow us 



N* N 




(a) (b) 

Figure 14: (a) pNN* vertex, (b) A''* — h bubble contributing to the selfenergy of the p meson, 
to calculate the selfenergy corresponding to the graph |l4) and we get 



In eq. (^), {IF) = 2 and {SF) = | are isospin and spin factors respectively, n{q) stands 



3 

for the occupation number of the nucleonic hole and £b{<1) is the energy of the baryon B 
involved, 

^2 



£B{q) = + g 2 ^Mb + (59) 



The integration can be performed analytically, and we are left with 



nf-^P) = -9hNpUN*{P) 



Un^{P) = , , (60) 

where we have defined the Lindhard function for the A^* — h excitation, Um*{P)- In both 
eqs. (|58|) and (|60| ) we have omitted the decay width in the N* propagator, although it is 
included in the calculation (see Appendix). 

Following the structure of the p meson propagator from the T matrix discussed in the 
previous section, now it is easy to include a selfenergy term coming from a pNN* coupling, 
just adding it to the one already existing from the coupling of the p meson to vrvr pairs, 
nucleons and deltas: 

+ (6i) 

5 Results and discussion 

First of all we show the dependence of our results on parameters that enter the calculation 
of the pion selfenergy, one of the essential ingredient in the renormalization of the p up to 



this point. In Fig. 15 we show the vrvr — > tttt scattering amplitude vs invariant mass. The 
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calculation has been performed neglecting the A^* — h excitation. The g' parameter in the 
pion selfenergy is set to 0.7. In the figure we plot the real and imaginary parts for the free 
case and at normal nuclear matter density. Compared to the free case, the imaginary part 
shows a clear broadening which is expected since there are new decay channels. At the 
same time one can observe a shift of the peak of the p distribution to higher energies. The 
resonance shape of the distribution remains at p = po and the zero of the real part of the 
amplitude also moves to higher energies. At normal nuclear density, the shift amounts to 
about 30 — 40 MeV and the width at half maximum is somewhat more than 200 MeV to 
be compared to the free value of 150 MeV. We have changed the irNN, ttNA form factor 
parameter A between 900 MeV and 1100 MeV. The differences found are small and this can 
give an idea of the uncertainties in the present results which we can expect from uncertainties 
in the pion selfenergy. Other uncertainties coming from details on the treatment of the A 
resonance will be discussed later. 




MI (GeV) 

Figure 15: Real and imaginary parts of the tttt —> nn scattering amplitude vs invariant mass. Solid lines 
stand for the free case. The others correspond to the calculation at p = po: dashed-dotted lines for A = 1 
GeV, long dashed and dotted lines for A — 0.9, 1.1 GeV respectively. 



We have also checked the dependence of the calculation on the g' Landau-Migdal pa- 
rameter appearing in eq. (p!6[). The results of this check are presented in Fig. 16, in which 
we show the real and imaginary parts of the scattering amplitude at p = pQ. The free case 
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is also displayed for reference. The effects of a varying g' are noticeable in the resonance 
region whereas the resuhs are insensitive for lower and higher energies. Neither the position 
of the resonance nor its width are much changed with variations of g' in the standard range 
of values [0.6,0.8]. We find a fluctuation of the peak position of around 3 % of the total 
mass at p = po- 




MI (GeV) 

Figure 16: The same as in Fig. ^ for p = po and several values of g' . Dashed-doted lines stand for 
g — 0.7, long dashed and dotted lines being for g' — 0.6, 0.8 respectively. 



The results on the mass shift presented so far are in qualitative agreement with other 
works which find a moderate shift of the p meson mass to higher energies at finite densities 
by using different approaches |jl6|, |l^] . 

The treatment of the A(1232) resonance is another issue with remarkable impact on 
the spectral behaviour of the p meson in nuclear matter. The works of refs. |16, where 



the main ingredient is the medium pion selfenergy mostly driven by the coupling to A — /i 
excitations find also a broadening of the p meson. We have found that the shape of the 
spectral distribution at low energies depends significantly on the A resonance properties 
specially its decay width. For instance, a simplified treatment in which a constant width 
is assumed gives a bump at around 500 MeV for p ^ po- This bump disappears when an 
energy dependent width is used (see also [^, Fig. 3.13). Nonetheless, in some calculations 
having explicit energy dependence on the A width the bump appearing at low energies in 
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the spectral function of the p meson, although weakened, still remains. 

Next we include in the calculation the A'^*(1520) contribution to the p meson selfenergy 
in nuclear matter, as described in section 4. Our model does not try to be complete since 
many other resonances should be included. A much detailed work along these lines can be 
found in ref. |^2|. Our aim is simply to estimate the effect of these type of decay channels 
on our previous results. The tttt amplitude is plotted in Fig. ^ for p = pQ with and without 
the N* — h contribution, together with the free case for reference. The possibility of exciting 
a N* — h pair is reflected in the structure arising in the amplitude around 550 MeV. As 
a consequence a sizeable amount of strength appears at energies below the p meson mass. 
Our results agree qualitatively with other works in which the p meson is allowed to couple 



to baryonic resonances [22, ^ 




MI (GeV) 

Figure 17: Real and imaginary parts of the tttt tttt scattering amplitude. Long dashed (dotted) lines 
correspond to the results without (with) the effects of the coupling pNN* at p = po, and solid lines stand 
for the free case. 

The same behaviour discussed above is reflected in the in-medium spectral function 



Sp = —^ImDp, which we show in Fig. 18 for the free case and several values of nuclear 
matter density. The first effect to be commented is that the peak of the p distribution 
broadens substantially as the density is increased, its width being about 200 MeV for p = pQ. 
This is accompanied by a clear shift of the position of the resonance to higher energies (in 
agreement with many other works, see Q for reference) that amounts, as reported before, 
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MI (GeV) 

Figure 18: Spectral function of the p meson for several values of nuclear matter density. Lines are as 
follows: solid, dashed and dotted stand respectively for p = 0, po/2 and pQ. 



to about 30 — 40 MeV for normal nuclear matter density. Together with the p peak we 
can observe a resonant structure corresponding to the excitation of the A^*(1520) appearing 
around 500-550 MeV. This increases noticeably the amount of strength towards low energies, 
which goes in the right direction in what concerns the description of dilepton radiation 
spectra in that energy region Q. The strength of the A^*(1520) peak increases with density, 
which is expected since this is an effect linear in density in our model. It is also noticeable 
that the position of this structure is slightly pushed downwards in energy with increasing 
density, as a consequence of the interference with the other contributions to the p selfenergy. 
In this respect we agree with refs. [24, 25] where it was pointed out that because of the 
coupling to baryon resonances, the in medium vector meson strength is split into a meson 
like mode, which is pushed up in energy and a resonance-hole like mode, which is pushed 
down in energy. 

The results would change if one takes into account that the properties of baryonic 
resonances depend themselves on the rho meson medium spectral function, as it is indeed the 
case for the A^*(1520) since it decays, for instance, into N pI'ktt]. Therefore a selfconsistent 
treatment is advisable. This has been done in |2^, where a melting of both A^* and p meson 
structures is found at high densities so that none of them are distinguishable any longer. As 
a consequence, much strength is spread to lower invariant masses. Note however that this 
result cannot be extrapolated directly to our case since we obtain a weaker contribution of 
the A^* — h excitation and our in- medium p is also narrower than in |22]. 

Our results reflect the most important characteristics of the medium corrections and 
these are shared qualitatively by many other works in which different approaches to the 



problem have been used ||l^, 16, 17, 18, 25 1. 



23 



6 Summary 



We have analyzed the problem of the p properties in a nuclear medium from a new per- 
spective, using a chiral unitary framework to describe the p meson in free space and in the 
medium. 

The approach starts from the lowest order chiral lagrangian involving the pseudoscalar 
mesons plus those lagrangians coupling vector mesons to pseudoscalar mesons. The unita- 
rization for meson- meson scattering is done following closely the idea of the N/D method 
in coupled channels and leads to an iteration of the tree level diagrams obtained from the 
chiral lagrangians plus the exchange of bare vector mesons. The free parameters of the 
theory are obtained by matching the results to those of chiral perturbation theory for the vr 
and K form factors plus the position of the p peak. Once this is done the theory reproduces 
the experimental form factors and the scattering amplitudes in the vector sector with high 
accuracy up to 1.2 GeV. 

The treatment of the p in the medium has led us to modify the previous approach, 
which relies upon the on-shell amplitudes provided by the lagrangians and dimensional 
regularization, in order to properly incorporate pion selfenergies and vertex corrections in 
a gauge invariant way. Whereas in vacuum it has been found that a formalism keeping 
tadpole diagrams and full off-shell dependence of the pvrvr vertex is equivalent to the one 
followed in ref. |32] where only the on shell part of the /jvrvr vertex is kept and no tadpole 
diagrams are considered, in the medium one should explicitly include tadpole terms in order 
to preserve gauge invariance. 

The inclusion of medium effects proceeds by dressing the pion propagators with a proper 
selfenergy and including vertex corrections. In this way the main p decay channel into vrvr 
is modified and new decay channels are opened. The main visible effect of the nuclear 
medium is a broadening of the p width as the nuclear density increases. Simultaneously we 
observe that there is a small shift of the peak to higher energies moderately dependent on 
the uncertainties of the pion selfenergy. 

We have also added another source of medium p renormalization including the coupling 
of the p to the A^*(1520) — h component. This part is not novel and has been studied 
elsewhere, but we include it for completeness in order to show how our previous results are 
affected by this channel. We observe that an important source of strength appears in the p 
spectral function at low energies and that the previous peak of the p distribution is slightly 
pushed to higher energies. Yet the global effect is still a shift of the p spectral function 
peak to higher energies of about 30 — 40 MeV at /j = po- Our results would thus agree with 
those where there is not much shift of the p peak but however produce an extra strength 
in the isovector- vector channel at lower energies than the free p mass, although there are 
quantitative differences. 
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Appendix: Formulae for the N* — h Lindhard function 



We give the analytic expression for the Lindhard function j|3^ corresponding to the A*"* — h 
excitation, quoted in eq. (60): 



3 Mn* 



2'^ Pkp 



z + -{l-z^) log — ^ + z^ + -{l- zl) log -^—^ 



(62) 



where kp \s the Fermi momentum, and here P stands for the modulus of P. Terms de- 
pending on z, Zc come from direct and crossed diagrams, respectively. These functions are 
defined as 



Pkp 
Pkp 



(63) 



with w/j = Mtv* — M]\f. Eq. (63) already includes the information of the A^* decay width, 
which can be found i n p^ . 

The limit of eq. (q^) when P ^ is given by 



UN*iP^,0;p)=p 



^ + 



1 



(64) 



and it is the expression that we use in our calculation as we stay in the CM frame of the p 
meson. 



Appendix 2: Other medium corrections not included in the 
text 

In section 3 we have described how the presence of nuclear matter modifies the propagation 
of the p meson by modifying the properties of the pionic cloud. In addition to the pion 
selfenergy, a set of diagrams involving the pirNN contact term were also considered to 
preserve gauge invariance (see Fig. ^). Minimal substitution, which provides the interaction 
lagrangians for the mesons and baryons involved [0, 18 1, gives rise to some other diagrams 



that we have not included in the text because (i) they are higher powers of density (in 
the calculation only corrections of order p are considered), (ii) they vanish for vanishing p 
momentum or because of symmetry reasons, or (iii) they are negligible compared to other 
corrections. Now we focus on some of them and outline the derivation of some expressions 
that lead to our conclusions. For this purpose we classify the graphs in two sets, depending 
on the number of pion propagators (zero or one pionic lines), as it is shown in Fig. |T9| . 

No pion lines 



The first vanishing contribution corresponds to the diagram labelled as |19|a, in which a p 
meson directly couples to & p — h excitation. This mechanism does not contribute for a p 
meson at rest, since then both the particle and the hole have the same momentum. The 
similar diagram with A — /i excitation involves a vanishing vertex in the non relativistic 
approximation, for a p at rest. 
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(a) 





Figure 19: Medium correction graphs not included in the calculation; sets (a) and (b) involving zero and 
one pion lines. pNN, pAA vertices are displayed as dots. 



One pion line 

Here we distribute the graphs in three subsets: both p mesons coupled to either a particle 
or a hole line (Fig. |l9|b.l), one p meson line attached to the particle line and the other one 
to the hole line (Fig. [l9|b.2) and both p mesons involved in a ppNN contact vertex (Fig. 
He). 

Let us calculate the contribution of diagrams [l9|b.l. to the p meson selfenergy. We focus 
on the case in which both incoming and outgoing p mesons are coupled to the particle line 
(in Fig. the momentum labelling is shown). Applying the Feynman rules derived from 



the interaction lagrangians of [18| we get the following expression: 



Gh{p)Glip + q)Gp{p' + - V^,p + q)^ (65) 

where Gp (Gh) stands for the propagator of a particle (hole) and K = {^Ydfii^-^Y ■ In 
obtaining the previous expression, the pNN vertices provide a {p + q)i{p + q)j factor. The 
terms linear and quadratic in p, when performing the integration over the momenta in the 

4 

fermionic loop give rise to contributions of order p s and higher powers in the nuclear density. 
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(p-'+q-.p + q) 



\(q'' q) 



(^,0) 




(p + q-{s ,p + q) 



(pVq^p + q) 



(fi,0) 



(pVq'',p + q) 



(pVqVjs ,p + q)/k 



(^,0) 



(q°. q) 




(pVq^p + q) 



(^,0) 



Figure 20: Momentum labelling of the diagram with two p mesons coupled to the particle line. 



We neglect those contributions and then eq. (65) follows. This result can be expressed in 
terms of the ordinary Lindhard functions. The product of fermion propagators in eq. ( |65[ ) 
can be rewritten as 



d ( 1 



da \ ^/s + a 

•,2 



d r -1 



+ (7° - a,p + q) - Gh{p)Gp{p'^ + g° + ^s,p + 

I „o 



da \ ^/s — a 



Gh{p)Gp{p^ + - a,p + q) - Gh{p)Gp{p^ + - V^,p + 



(66) 



where the limit a — > is assumed. Commutation of this limit and the da derivative with 
the integration over the fermionic loop momenta allows us to write 



(67) 



where the U function, the Lindhard function for a forward propagating p — /i bubble, stands 
for 



UN{q) = -2i 



d'^p 



Gh{p)Gp{p + q) 



and satisfies 



(2vr)4 

2[UN{q^q)+UN{-q'',q)] = UN{q\q). 



(68) 



(69) 
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We have compared numerically the size of this contribution to any of the ones included in 
Fig. |6| and we find that it is around a hundred times smaller in the range of invariant mass 
under study. Including this contribution in the calculation has no visible effect and therefore 
we neglect it, and so we do with the terms containing A baryons instead of nucleons. 

In this derivation we have neglected all terms linear or quadratic in the hole momentum, 
since they generate contributions of higher order in density. Following this criterium, we 
can discard other contributions from diagrams with a p meson coupled to a hole line. In 
particular this is the case of the partners of the graphs shown of Fig. with both p mesons 
coupled to the hole line. The same happens to the diagrams of the type of the one shown 
in Fig. [lib. 2. 



We are left with the evaluation of the ppinT vertex correction shown in Fig. 19:, where 
although not displayed, the ppNN vertex can also be inserted in the hole line. The mo- 
mentum labelling follows the one in Fig. |2^. Using the Feynman rule for the ppNN vertex 
derived in section 3.6.3 we have 



(70) 



where here G{p) denotes the whole nucleon propagator with both the particle and the hole 
pieces. It is possible to rewrite G{p + q)^ as 

G(p + g)2 = a„G(p° + gO-a,p-*+gl, (71) 

with a — > 0. The integration of the fermionic loop is then given by 

I = I ^G{p)G{p + qf = d.j ^^G{p)G{p' + q'- a,p + q) = 

'-daU{q'> -a,q)\a=o, (72) 

which is an odd function of q^ since the Lindhard function is even in q^. As a consequence 
the dq^ integration vanishes and the contribution of these terms is null. Similar arguments 
for the case of a A in the particle line also lead to a null contribution. 
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